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The highly sensitive method to search for WIMPs dark matter particles is proposed. An 
array of thin Nal(Tl) plate has the great selectivity for distinguishing the WIMPs events 
and background ones. The principle of signal selection for WIMPs is described. The high 
sensitivity for SD (spin-dependent) type WIMPs is expected by applying multi-layer system 



(N 
> 

Q\ ' of Nal(Tl) detector 

(N ' 

CO ■ KEYWORDS: Thin scintillator, WIMPs search, Dark matter 

o 
in 
o 



1. Introduction 

A method of highly sensitive measurement is strongly needed to study nuclear and particle 
rare processes. The need for highly sensitive radiation detector has been increased in the 



fields of particle astrophysics and nuclear physics. Particularly, the search for dark matter 
candidates and the search for neutrino-less double beta decay need to select the signal events 
from the huge number of background ones. In the case of WIMPs search, the groups such as 
ELEGANT, 1 DAMA 2 applied the huge volume of Nal(Tl) scintillator. The new project with 
the huge volume scintillator are proposed with hundreds kg of scintillators by XMASS 3 group, 
ZEPLIN group 4 and so on. 

The many groups are trying to enhance the sensitivity by developing a large volume and 
low background detectors. However, to reduce the background events becomes more difficult 
for the larger detectors because of its poor position resolution. We propose highly selective 
detector system to extract the signals of WIMPs. The segmentation of a detector results good 
position resolution, which enhances the sensitivity for WIMPs. 

The three types of interaction between WIMPs and nucleus has been proposed by many 
theorists. One is spin-independent (SI), the second is spin-dependent (SD), and the third 
one is spin-dependent inelastic scattering (EX). The SI interaction cross section depends on 
the square of nuclear mass number A 2 . On the other hand, the SD one depends on the 
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nuclear spin-matrix element A 2 J(J + 1). The EX interaction cross section is strongly related 
to the spin-dependent (SD) interaction cross section. Consequently, the event rate of EX and 
SD are related each other, moreover, the relationship is enhanced by measuring the annual 
modulation. The highly selective detector which measures both EX and elastic scattering 
contributes much for the WIMPs search. 

It consists of an array of thin Nal(Tl) scintillator plate. The event selection power is 
largely enhanced by applying the compact modular detector system. The size of each modules 
is designed in order to enhance the discrimination between the signal and the background. 

2. Signal Selection by Space and Time Correlation (SSSTC) 

The WIMPs signal is expected to make the localized event both in space and in time. On 
the other hand, the background events make the diffused and correlated ones. Consequently, 
the spatial correlation and timing correlation is quite important information to extract the 
signal event from the large number of background events. Fig.l(A) shows the schematic spatial- 
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Fig. 1. Schematic drawing of the principle of SSSTC. A detector is spatially segmented by the vertical 
dashed lines. 



timing correlation of the event in the case of WIMPs-nuclear inelastic excitation. The number 
1,2,3 and 4 stand for the order of events occurred. The event indicated "1" is the signal event 
for WIMPs-nucleus interaction, and the ones "2" and "3" are background events. 
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In the case of elastic scattering, the events cannot be the spatially diffused event, and 
must be the timing isolated event because the range of recoil nucleus (only "R" in Fig.l) is 
much shorter than the thickness of each module. 

In the case of inelastic excitation of 127 I nucleus, the 57.6keV 7 ray is simultaneously 
emitted. Such the signal events for WIMPs- nucleus interaction produce the characteristic 
events. The recoil nucleus "R" produces a continuous energy spectrum in the low energy 
region and the low energy 7 ray "7" in Fig.l is observed in the another modules. 

Fig. 1(B) shows that the RI background event which makes a similar event. The (3 — 7 
successive transition makes similar event to the signal ones. The most serious origin of such 
a fake event is due to 210 Pb. The radioactivity of 210 Pb in Nal(Tl) was about lOmBq/kg, 
which is one thousand times larger than that of 226 Ra. 5 The background events due to 210 Pb 
is tagged by the following j3 decay event by 210 Bi and eliminated. The accidental background 
events suffers the event selection in the case of a large volume detector. Many groups has 
tried to achieve good position resolution by means of pulse shape discrimination or timing 
discrimination or pulse height dependence. However, it is quite clear that the segmentation 
of the detector is achieved most effectively by means of small modules of detector. In the 
following sections, the effectiveness of event selection of WIMPs- nucleus inelastic excitation 
(section 3) and the reduction power of 210 Pb background (section 4) by means of segmentation. 

3. Event selection by means of segmented Nal(Tl) 

In this section, the selectivity of EX events by highly segmented scintillator is discussed. 
Since the attenuation coefficient of 57.6keV 7 ray is as large as 6.83cm 2 /g, each detector 
module must be sufficiently small. We determined the design of the Nal(Tl) scintillator so as 
to detect the 7 ray by the next module with enlarging the volume of the detector. The Nal(Tl) 
crystal with 5cm x 5cm in area and less than 2mm in thickness was considered by performing 
Monte Carlo simulation. 

The conditions which has been considered in the simulation are listed below. 

(1) The thickness of Nal(Tl) crystal was varied 0.1mm to 1.5mm. 

(2) The wider area of the Nal(Tl) crystal was covered with ESR™ (Enhanced Specular 
Reflector)reflector sheets. 

(3) Four light guides made of acrylic plate were placed on the thinner sides of Nal(Tl) crystal. 

(4) A reinforcement was glued on ESR™ sheet. 

(5) Ten modules were piled up. 

(6) The inelastic excitation occurred in the 5th module. 

(7) The simulation was performed with 1 million events, for EX and 210 Pb, 3 million events 
for 210 Pb. 
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(8) The energy resolution of the thin Nal(Tl) was assumed as large as 19% in FWHM (Full 
Width Half Maximum) at 60keV. 

The drawing of one module is shown in Fig.2. Optical cement is glued between Nal(Tl) and 
optical window. 
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Fig. 2. The schematic drawing of a Nal(Tl) module. 



The analysis was performed in order to extract the EX events. The number of the events 
in the energy range between 47keV and 67keV (~57.6keV±10keV) which were detected by 
4th or 6th module was divided by the number of total simulated events (1 million events) 
to estimate the 'coincidence efficiency'. Fig. 3 shows the thickness dependence of coincidence 
efficiency. 

The coincidence efficiency has a maximum between 0.3mm and 0.5mm. In the case of 
the Nal(Tl) is thinner than 0.3mm, the 7 ray goes through the next module, which makes a 
spatially diffused event. On the other hand, the coincidence efficiency does not decrease rapidly 
when the Nal(Tl) becomes thicker. However, since the thicker module makes larger event rates 
of background events, the background events suffers the sensitivity. The background reduction 
power will be discussed in the next section. 

Recently, the suitable thin Nal(Tl) crystal with the wide area in collaboration with Horiba 
Ltd. The thin Nal(Tl) detector whose dimension was 5cm x 5cm x 0.05cm was successfully de- 
veloped. 6 The thin Nal(Tl) plate with 0.05cm in thickness is suitable for SSSTC measurement 
of WIMPs search. The performance of the thin Nal(Tl) was tested by low energy photons. 
The good energy resolution of 19% at 60keV and low energy threshold of about 2keV was 
obtained by thin Nal(Tl) detector. 6 
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Fig. 3. The thickness dependence of coincidence efficiency. 



4. Background reduction by means of segmented Nal(Tl) 

The background events for WIMPs search are mainly due to the internal RI's (radioactive 
isotopes) which are contained in the sensitive volume of the detector. The external RI's which 
are contained in the surrounding materials also make serious background. Since the back- 
ground events due to external RI's are mainly due to Compton scattering of high energy 7 
rays, the background events are effectively reduced by the highly segmented detector system. 
Moreover, the internal RI's are much more serious than external ones because of the large 
geometrical efficiency. Consequently, the background events due to only the internal RI's were 
considered in the present estimation. 

The most serious origin of background is 210 Pb. 210 Pb emits low energy (3 ray (E max = 
16.5keV) and low energy 7 ray (Ej = 46.5keV). The simultaneous emission of the low energy 
(3 ray and the 7 ray makes a similar event to a signal event of EX. The branching ratio that 
210 Pb decays to the first excited state is 84%. The transition from the first excited state to the 
ground state is mainly by internal conversion, consequently, the probability of 7 ray emission 
is 4.25%. 

The normal Nal(Tl) crystal contains a large amount of 210 Pb. For example, the concen- 
tration of 210 Pb was 8.4mBq/kg in the Nal(Tl) of ELEGANT V, while the one of 214 Pb which 
is the parents of 210 Pb was as small as a few tens of /xBq/kg. The radiation equilibrium in 
U-chain in Nal(Tl) is often broken. 7 ' 8 The progeny after 222 Rn in U-chain were concentrated 
from the air, water and raw materials. 

It is quite difficult to reduce the concentration of 210 Pb, however, its contribution to 
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background is efficiently reduced by SSSTC analysis. In the segmented detector module, the 
progeny of 210 Pb decays after its decay and detected in the same module. Suppose the 1kg 
detector is divided into n modules which contains &Bq/kg of 210 Pb. The decay rate in one 
module is simply expressed as 6/nsec _1 module _1 . The sequential events which occur during 
the interval AT are rejected as the decay chain of 210 Pb^ 210 Bi^ 210 Po. In this case, the 
accidental decay due to another 210 Pb must not occur, thus the mean time interval due to 
210 Pb must be longer than the time interval. Thus, the condition 



is needed. The reduction power for BG events is enhanced with setting a longer time interval 
AT. The time interval AT = N x T 1/2 , with N = 2 ~ 3 is adequate for SSSTC analysis. 

The present WIMPs search project aims to develop the highly pure Nal(Tl) crystal whose 
concentration of 210 Pb is as small as O.lmBq/kg. The target concentration of 210 Pb is ade- 
quately realistic purity. In this case, segmentation n is calculated as 



The proposed Nal(Tl) module whose dimension of 5cmx 5cm x 0.05cm corresponds to 217 
segmentation. The high sensitive measurement is achieved by the highly segmented detector 
system with the feasible purity of background radioactivity. 

5. Estimation of the sensitivity for inelastic excitation 

The expected sensitivity for inelastic scattering of WIMPs- 127 I will be discussed in this 
section. The expected background events were estimated by means of Monte Carlo simula- 
tion. The serious background origins, 214 Pb, 214 Bi, 210 Pb and 210 Bi in were considered. The 
contribution of 40 K is expected to be small because of extremely small radioactivity (less than 
0.2ppm in nat K) in the present Nal(Tl). The simulated concentration of 214 Pb and 214 Bi was 
O.OlmBq/kg which was precisely measured by delayed coincidence method. 

While the simulated concentration of 210 Pb was O.lmBq/kg, which was 1/100 smaller than 
the present concentration. The radioactive equilibrium in U-chain often breaks because of the 
concentration of 222 Rn in the air. Our previous work found that the concentration of 210 Pb is 
almost three orders of magnitude larger than the one of 214 Pb and 214 Bi. The contamination 
of 210 Pb is reduced by controlling the air during the creation processes of Nal(Tl) crystal. 
The Monte Carlo simulation was performed under the condition that O.lmBq/kg in Nal(Tl) 
crystal. 

The SSSTC analysis was performed for the simulated background data. A event in the 
energy region 57.6keV±10keV was selected for the inelastic excitation of 127 I. After the anal- 
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ysis, the accidental coincidence events were reduced by setting the timing correlation. The 
successive events during the specific time were concluded to be the decays belonging to the 
decay chain. In the case of the chain between 214 Pb and 214 Bi, the event which occurs 1 second 
and 3600 seconds after the preceding event is removed as the background one. The fraction 
of removed events is estimated as 

p3600 



r© 



dt = 0.997. 



In the case of the decay chain from 210 Pb to 210 Po, the 82.3% of the background events are 
removed with setting the time difference being 12.5days (2.5^/2). 

The estimated background energy spectrum is shown in Fig. 4. The energy resolution 
which was measured by actual thin Nal(Tl) was considered in this simulation. The events 
between 2keV and lOkeV was summed to estimate the sensitivity for WIMPs-nucleus inelastic 
excitation. The energy spectrum after performing the SSSTC analysis has the big contribution 
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Fig. 4. The expected background energy spectrum. The closed circles and open circles mean the 
singles event rate and the event rate after performing the SSSTC analysis. The analysis process 
is described in text. 



by 210 Pb. A big bump below 20keV of the open circle is due to the beta ray from °Pb 
(E max = 16.5keV). 

The expected event rate due to the radioactive contamination which is contained in the 
Nal(Tl) crystal was calculated as 1.7x 10~ 2 /kg/day. The assumed energy threshold was 2keV 
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in electron equivalent and the event rate between 2keV and lOkeV was integrated. The energy 
threshold, 2keV, was already achieved by thin Nal(Tl) scintillator. 6 The statistical accuracy 
is improved as the increase of the modules. The expected background rate and the upper limit 
(90% C.L.) of the event rate are listed in table I. 



Table I. The expected event rate for multi layer detector system of Nal(Tl) plate. The upper limit 
on the event rate at 90% C.L. are shown. 



Number of modules 


Total mass (kg) 


BG rate(kg- 1 y~ 1 ) 


Upper limit (ifemkg'May" 1 ) 


16 


0.0734 


0.45 ± 0.67 


5.2 x 10" 2 


256 


1.17 


7.2 ±2.7 


8.1 x 10~ 3 


1024 


4.70 


28.9 ±5.4 


4.0 x 10~ 3 


2176 


9.98 


61.3 ±7.8 


2.8 x 10~ 3 



The expected sensitivity for the cross section of inelastic excitation is calculated easily by 
the formula, 

VEX,lim = — 2 X R Hm- (3) 

NtPo(v) f\F(q)\ e 

Where m x and are the mass of WIMPs and the number density of target nucleus ( 127 I), 
po = 0.3GeV/cm 3 is the local halo density, (v) = 230km/sec is the mean velocity of WIMPs. 
The phase space factor / = Pf/pi was calculated numerically. The form factor F(q) was used 
for the spin-dependent form factor. 9 The coincidence efficiency e = 0.21 is discussed in the 
previous section. 

The other background from out of the detector should be considered. The background is 
mainly due to high energy gamma rays from photomultiplier tubes (PMT). The high energy 
gamma rays are reduced because they interact in many modules through multiple Compton 
scattering. Since the multiple Compton scattering events are effectively reduced, the BG's due 
to external origins of the detector is ignored in the present report. 

The cross sections of EX and SD are closely related by the nuclear spin matrix element. 
The case of elastic scattering, the nuclear spin-matrix element, X 2 J(J + 1) has the large 
model dependence. 9 The spin-matrix element using the calculation by means of single particle 
shell model was applied to our estimation. 9 However, the spin-matrix element for 127 I cancels 
out when one renormalizes the sensitivity to WIMPs-proton cross section. The spin-matrix 
element for proton has no model dependence and the value is A 2 J( J + 1) = 0.75. 

The case of inelastic scattering the nuclear spin-matrix element, ^fj+i ~ff 41 * s ex P er i~ 
mentally deduced from the nuclear transition probability, where J and J* are the total angular 
momentum of the ground state and excited state, and [x p is the magnetic dipole moment of 
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proton. Consequently the precise exclusion plot with small model dependence is obtained. 
The sensitivity of the cross section au mtP - x is calculated simply, 

^ m , p _ x =[A 2 J(J + l)] p ^|f^±^) 
L Jp m z N \rriN + m x J 

In the present estimation, the parameters in the formulae are listed in table II. 



M M i 



2 J* + 1 
2J + 1 



<riim,EX- 



(4) 



Table II. The parameters which are used for the present estimation. 



Parameter 


Value 


Ref. 


m p 


938GeV 


10 


m N 


118.18GeV 


10 


J 


5/2 


10 


J* 


7/2 


10 


\M M i\ 2 


0.10 


11 




2.79 


10 



The estimated sensitivity for SD type WIMPs is shown in Fig. 5. The high sensitivity will 
be expected by the small amount of Nal(Tl) crystal. 

6. Concluding remarks and future prospect 

The WIMPs search by means of highly segmented Nal(Tl) scintillator was proposed in 
the present work. It is important to investigate both the elastic scattering and the inelastic 
scattering of WIMPs-nucleus interaction. The ratio between elastic scattering and inelastic 
scattering provides the quite important information for WIMPs nature. We showed the thin 
and wide area Nal(Tl) plate was suitable for study the WIMPs interaction. 

Since Nal(Tl) scintillators have been well developed, the production of highly pure and 
large volume detectors are provided easily and with low cost. The thin Nal(Tl) plate has been 
developed for this project by Horiba Ltd. successfully with low cost and with high purity. The 
Nal(Tl) array needs the large number of channels for data acquisition (DAQ) system. Recently, 
DAQ systems for the large number of module (1024-2048 channels) are extensively developed 
in the fields of high energy experiment. The low noise measurement should be performed for 
the present work. Since the developed thin Nal(Tl) has the low energy threshold (~2keV), 
the noise from PMT will not be a big problem. However, the pulse shape analysis will be 
performed in order to monitor the PMT noises. 

We showed that the thin Nal(Tl) array has excellent sensitivity for WIMP-nucleus inelastic 
excitation by only 10kg of Nal(Tl). The highly segmented detector system has the great 
advantage to SD-type WIMPs candidates. 
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Fig. 5. Thick lines are the expected sensitivity for SD type WIMPs. The solid, long-dashed, sort- 
dashed and dotted lines are the expected sensitivity by 16, 256, 1024 and 2176 modules of Nal(Tl) 
array. Thin lines arc results of CRESST (solid line), DAMA LXc (long-dashed line), ELEGANT 
V at OTO (short-dashed line) and the expected sensitivity by NAIAD (dotted line). 
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